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Abstract 

We present observations of the first ten degrees of longitude in the Mopra carbon monoxide (CO) survey 
of the southern Galactic plane ( Burton et al.||2013 ), covering Galactic longitude I = 320-330° and latitude 
b = ±0.5°, and I = 327-330°, b = ±0.5-1.0°. These data have been taken at 35 arcsec spatial resolution 
and O.lkrns”^ spectral resolution, providing an unprecedented view of the molecular clouds and gas of the 
southern Galactic plane in the 109-115 GHz J = 1-0 transitions of ^^CO, ^^CO, and C^^O. Together 

with information about the noise statistics from the Mopra telescope, these data can be retrieved from the 
Mopra CO website and the CSIRO-ATNF data archive. 


Keywords: Galaxy: kinematics and dynamics - Galaxy: structure - ISM: clouds - ISM: molecules - radio 
lines: ISM - surveys 


1 INTRODUCTION 

Molecular clouds comprise the densest regions of the in¬ 
terstellar medium, and are the nurseries in which stars 
form. In order to determine how molecular clouds them¬ 
selves form, surveys of the three principal forms of car¬ 
bon (i.e. C+, C and CO) are needed to trace the gas in 
its ionised, atomic and molecular forms respectively. In 
particular, surveys along the Galactic plane, where the 
great majority of molecular clouds lie, provide us with 
a picture of the molecular structure of our Galaxy. 

The Mopra southern Galactic plane CO Survey is de¬ 
signed to map the distribution and dynamics of the car¬ 
bon monoxide (CO) molecule along the Fourth Quad¬ 
rant of the Galaxy, from I = 270-360° with b = ±0.5° 
( Burton et al.|[^13[ Paper I). Using the Mopra radio 
telescope, the J = 1-0 line of CO at 2.6 mm is be¬ 
ing mapped at 35 arcsec spatial and 0.1kms“^ spectral 
resolution, in order to determine where giant molecu¬ 
lar clouds (GMCs) lie in the Galaxy and to explore 


the connection between molecular clouds and the ‘miss¬ 
ing’ gas inferred from gamma-ray observations. The CO 
data will be combined with similar survey data col¬ 
lected of the sub-millimetre neutral carbon line (C) 
taken with the Nanten2 telescope in Chile, and the tera¬ 
hertz neutral and ionized carbon lines (C and C+) from 
the balloon-borne Stratospheric Terahertz Observatory 
2 (ST02) and the ice-bound High Elevation Antarctic 
Telescope (HEAT), both in Antarctica. With these sur¬ 
veys it will be possible to determine the distribution 
and motion of the principal forms of carbon in the in¬ 
terstellar medium. 

Data Release 1 (hereafter DRl) covers the first ten 
square degrees of the survey, from I = 320-330° with 
b = ±0.5°, along with an additional three half-square 
degrees over I = 327-330°, b = ±0.5-1.0° (see Eigures[^ 
&|^. The G323 sightline has previously been published 
in Paper I, however, it is discussed here because the 
G323 data was taken in the summer months and com¬ 
pares less favourably to the other data cubes, which 
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Spitzer 8 micron, 24 micron, and Mopra data 
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Figure 1. Three-colour image of the Galactic plane from the Spitzer GLIMPSE 8 p,m (red; [Ghurchwell et al.||2009[ [Benjamin et al.| 
|20Q3| > and MIPSGAL 24/rm (green; [Carey et al.|2009| surveys, with a moment 0 map calculated over v = —120 to 0 km s ^ from the 
Mopra ^^GO J = 1-0 survey (blue). The region enclosed by the white box, containing Galactic longitudes I = 320-330° and latitudes 
b = ±0.5°, as well as I = 327-330°, b = ±0.5-1.0°, is that published here as Data Release 1 (DRl). 
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Figure 2. The moment 0 map calculated over v = —120 to Okms ^ from the Mopra ^^GO J = 1-0 Southern Galactic Plane survey. 
The region illustrated is that published here as Data Release 1 (DRl). 


were obtained in the Austral winter. The improved at¬ 
mospheric conditions make it possible to observe all four 
of the major CO isotopologues (^^CO, ^^CO, and 

C^^O), although the line is very weak and only 

detectable at a 2cr level at best, as measuring it typically 
requires longer integration times than our fast-mapping 
survey mode allows (see Paper I). The isotopologue ra¬ 
tios allow examination of the distribution of columns 
of molecular gas through GMCs ( |Heyer et al.||2QQ9 ), so 
that variance in the Xco factor, used to link the CO 
intensity to the H 2 column density in order to obtain 
molecular cloud masses, can be determined (see e.g 
latto, Wolfire fc Leroy||2Q13 ). 


Bo- 


In Section 2 of this paper (Paper II) the data are de¬ 
scribed, along with the various metrics used to quantify 
the noise and data quality. Section 3 will present some 
sample science results, however detailed science inves¬ 


tigations are left for focussed studies, such as Burton 


et al. (2014), which described a filamentary molecular 


cloud in the G328 sightline. The data are being made 
available at the survey websit^ and in the CSIRO- 
ATNF data archiv^ as they are published. 


^http://www.phys.unsw.edu.au/mopraco/ 
^http: / / atoa.atnf.csiro.au 
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2 DATA RELEASE 1 


The data presented in this paper were taken in March 
2011, and over the Austral winters of 2011-2013. 
The observations performed in 2011 (/ = 323-330°, 
6 = ±0.5°) used the UNSW Digital Filter Bank (the 
UNSW-MOPS) in its ‘zoom’ mode, with 4 x 137.5 MHz 
dual-polarization bands; in 2012 this was increased to 
8 bands, so that the data from I = 320-323°, b = ±0.5° 
and I = 327-330°, b = ±0.5-1.0° cover a larger spectral 
region. Very little emission is observed in these extended 
bands away from the Galactic centre, but it is possible 
that high velocity clouds could be detected. The full 
line parameters are given in Table 1 of Paper I; the 
DRl dataset has a spectral resolution of ^ 0.1kms“^ 
over at least 4096 channels for each of the four GO iso¬ 
topologues. 

The angular resolution of the Mopra beam at 


115 GHz is 33arcsec FWHM (Ladd et al. 2005); after 


the median filter convolution is applied during the data 
reduction process, the beam size is around 35 arcsec in 
the final data set. The extended beam efficiency, used 
to convert brightness temperatures into line fluxes, is 
^XB = 0-55 at 115 GHz (rather than the main beam ef¬ 
ficiency ?7 mb = 0.42; [Ladd et al.|[2005 ). 

There are a number of metrics used to assess the qual¬ 
ity of DRl. Figureshows histograms of the noise per 
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Figure 4. Maps showing the noise level (in [K] units) for and ^^CO, determined from the standard deviation of the continuum 

channels between 0 and +90kms“^ for each pixel. Note that the intensity scales differ, as the observations have ~ 3x higher 

noise than ^^CO. The higher noise in the G323 region (Paper I) is also evident. 
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Figure 3. Probability distribution of the noise level in the 
(left) and data cubes, CTcont, as determined from the stan¬ 

dard deviation in the continuum channels (in [K] units) for 
each pixel. The dotted distribution in each is the noise level in 
the I = 323-324° field (described in Paper I), highlighting the 
effect of observing in the summer months. 


channel, (Jcont, as well as the previously-published dis¬ 
tribution from the G323 sightline for both the 
and lines (the and lines have similar 

noise characteristics to ^^CO). These are calculated by 
finding the standard deviation of the channels outside 
of the range where line emission occurs, where Mlsr < 
— 120kms“^ and Mlsr > +20kms“^ (Figure 11 shows 
that some emission exists where Vlsr > +20kms“^, 
but it occurs in very few pixels and does not affect 
these statistics). Although the G323 data were taken 
in the summer months, the noise distribution is not sig¬ 
nificantly different from those data obtained in winter. 
While the tails of the histograms are ^ 0.5 K higher in 
the summer observations due to the poorer conditions. 


the mode values for and ^^GO are 1.5 and 0.7K 

per O.lkrns”^ velocity channel, in comparison to 1.3 
and 0.5K per 0.1kms“^ channel in winter. The 
line has higher noise than the ^^GO line (and those of 
the other isotopologues not shown here), as there is a 
molecular oxygen absorption line near 115 GHz, making 
the atmosphere inherently worse than it is for the other 
isotopologues at 109-112 GHz. 

Figure 1^ shows the la noise maps for the ^^GO and 
data cubes, determined from the standard de¬ 
viation of the continuum channels between = 0 and 
V = +90kms“^ in each pixel. The wider visible strip¬ 
ing in these maps is an effect of the scans in the I and b 
directions being taken with a variety of different system 
temperatures and observing conditions. Thinner stripes 
or spots of lower apparent noise occur where ‘bad’ pix¬ 
els, rows or columns were identified and interpolated 
over (the lower noise is a result of interpolating over 
several nearby pixels; see Paper I). As in Figurethe 
higher noise occurring in the G323 field is a result of 
poorer conditions in the summer months. 

The system temperature, Tgys, measures the level 
of the received signal from the source, sky, telescope 
and instrument. It is calibrated every 30 minutes with 
reference to an ambient temperature paddle that is 
placed in front of the beam. Figure shows histograms 
of the probability distribution of Tgys for ^^GO and 
^^GO, alongside their counterparts for the G323 sight- 
line. Again, the line has Tgys about a factor of 

two higher than the values for the other lines due to 
the inferior conditions. Maps of the system tempera¬ 
ture for all four of the GO isotopologues are shown in 
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Figure 6. Tgys images for (from top): ^^CO, and (which share the same 2 GHz band of the correlator), and C^^O, in 

units of (K) (as indicated by the scale bars). The striping pattern is inherent to the data set, resulting from scanning in the I and 
b directions in variable observing conditions. The darkest square degree region, I = 323-324°, was observed in March 2011 when the 
summer conditions were less favourable for CO observations. 
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Figure 5. Probability distribution of the system temperature, 
Tgys (in [K] units) for and ^^CO, in the data for each 

pixel, determined from the ambient temperature load paddle mea¬ 
surements. The dotted distribution in each is the I = 323-324° 
field described in Paper 1; the higher Tsys values highlight the 
effect of observing in the summer months. 


Figure]^ and share the same Tgys map as 

these lines are in the same 2 GHz band of the correla¬ 
tor. As before, the striping is a result of combining data 
from scans in both the I and b directions (minimizing 
artefacts from poor sky conditions), while spots occur 
where the data cube has been thresholded to remove 


that data with excessive Tgys values (see Paper I for 
more details). Scans with particularly high Tgyg values 
were repeated (these are also visible in Figure Has wide 
stripes of lower noise). 

Several of the data cubes were affected by contam¬ 
ination from faint molecular clouds in the sky refer¬ 
ence beam (of typical intensity ^ —0.5K in and 

^ —O.lKin^^CO; the other lines are generally too weak 
for sky emission to be visible). While tests of the sky ref¬ 
erence positions were carried out prior to survey obser¬ 
vations, these were quick observations and the faint sky 
signal was typically lost in the noise. In DRl this con¬ 
tamination has been removed by calculating the average 
spectra in a region of ^ 100 pixels where the sky con¬ 
tamination is uniform, such that there is no true emis¬ 
sion within ±3kms“^ of the negative emission profile. 
A Gaussian template is fit to this contaminating profile 
and then subtracted from the individual spectrum at 
each pixel in the full square degree. The central veloc¬ 
ity of the sky emission varies across each square degree, 
as the reduction process (see Paper I) does not control 
for any variation relative to the local standard of rest 
in the sky spectrum; this is compensated for by linearly 
varying the central velocity of the template spectrum in 
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Table 1 Sky reference beam positions for each square degree of DRl; the additional half-square degree regions from I = 327-330°, 
b = +0.5-1.0° use the same sky reference positions as their main survey counterparts. Several of the cubes were affected by contamination 
due to line emission at the sky reference position; the central velocity and intensity (in [K] units) of this contamination is 

also noted. 


Cube 

Reference 

Longitude 

Reference 

Latitude 

Contamination 
Velocity (kms“^) 

Contamination Peak 
Intensity (T^ [K]) 

G320 

320.375 

-3.000 



G321 

321.500 

-2.500 

-1.68 

-3.18 

G322 

322.500 

-2.000 

-56.11, -12.96 

-1.25, -0.39 

G323 

323.500 

-2.000 



G324 

324.500 

-2.000 

+0.90 

-0.78 

G325 

325.500 

-2.200 

-46.62 

-0.49 

G326 

326.500 

-2.000 

-19.31 

-0.55 

G327 

327.500 

-2.200 



G328 

328.500 

-2.200 

+2.68 

-0.42 

G329 

329.500 

-2.200 




both longitude and latitude (typically by < 1 km in 
either direction). Table gives the position of the sky 
reference beam for each square degree, along with the 
intensity and central velocity of the contamination in 
the average spectrum for the square degree. 

DRl has been made available at the survey website 
and in the CSIRO-ATNF data archive as a series of 
fits files for each isotopologue, each containing a square 
degree of the survey, centred at each half degree along 
the Galactic plane. Although the full range of veloc¬ 
ity space available for is [—1000,+1000] kms“^ 

(with file sizes ^ 800MB / square degree), it is ex¬ 
pected that most users would prefer to download the 
smaller (^ 300 MB / square degree) files covering the 
velocity range [—150, +50] kms“^ illustrated in this pa¬ 
per. No emission has yet been detected outside this 
range, although closer to the Galactic centre a wider 
velocity range will be included within the smaller files. 
Also available are the Tgyg, beam coverage and rms noise 
images for each line. The 10° x 1.5° cube used for the 
analysis in this paper was created using the imcomb task 
in miriac0 ( jSault, Teuben fc Wright|[l995 ). 


3 RESULTS 


Figures and show the average ^^GO and ^^GO 
line profiles for each square degree cube in DRl. The 
^^GO profiles from the lower spa tial (8^) and spectral 
(1.3kms“^) resolution survey of Dame, Hartmann & 


Thaddeus| ( |2001 ) are overlaid on the plots, showing that 


the structure of the profiles match well, although the 
intensity of the Mopra line peaks is ^ 1.4 times higher 
than the corresponding [Dame, Hartmann fc: Thaddens] 
(2001) intensities (see Table [^. The intensity ratio was 
determined using the average spectra for each square 


^ http: / / WWW. atnf.csiro. au / computing/software / miriad / 
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Table 2 The median and mode of the ratio of the Mopra CO 
intensities to those of the [Dame, Hartmann &: Thaddeus| ( |2001| ) 
survey, calculated using the average spectra in each square degree 
over the velocity range v = [—120, 20] kms“^. 


Cube 

Median 

Mode 

G320 

1.42 

1.40 

G321 

1.31 

1.27 

G322 

1.28 

1.28 

G323 

1.42 

1.37 

G324 

1.62 

1.61 

G325 

1.41 

1.37 

G326 

1.34 

1.31 

G327 

1.38 

1.27 

G328 

1.64 

1.60 

G329 

1.41 

1.37 

Average 

1.42 

1.39 


degree over the velocity range i; = [—120, 20] kms“^ 
where the majority of the emission lies, and has been 
found to be relatively consistent across DRl. [Wong et 


al. (2011) found the average of the Mopra extended 


beam efficiency (t/xb) had a typical rms deviation of 
^10% over their eight observing seasons (each about 
a month in duration), with much of the variation re¬ 
flecting changes in the system temperature scale due to 
instrument modifications. Their beam efficiencies were 
lower than those used to calculate Tmb in this paper; 
using these would only increase the ratio of the Mopra 
CO fluxes to those of [Dame, Hartmann fc Thaddensj 
( |2QQT] ). 

The integrated intensity maps in Figures [9}|^ high¬ 
light molecular clouds in the Mopra and ^^CO 

data cubes over each 10kms“^ interval between —110 
and +40kms“^ in units of Kkms“^, corrected for the 
beam efficiency tjxb = 0.55. These have been calcu- 
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Mean Profile - G320 



Mean Profile - G321 



Mean Profile - G322 



Mean Profile - G323 



Mean Profile - G324 



Velocity Vlsr (km/s) 

Figure 7. The CO line profiles, averaged over each square degree from I = 320-325°, in units of Tmb (K) (i-e. divided by the telescope 
efficiency, ryxB = 0.55). The solid line is ^^CO; the dotted line is (which has been multiplied by 3 and offset by —0.5K for 

clarity); these spectra are binned in velocity by 2 pixels (0.2kms“^). The dot-dashed line is the equivalent spectrum from the 

[Dame, Hartmann &: Thaddeus] ( |2001| > survey, and the vertical dashed lines are the velocities from which contaminating sky emission 
was removed (see Table [fj. 
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Mean Profile - G325 



Mean Profile - G326 



Mean Profile - G327 



Mean Profile - G328 



Mean Profile - G329 



Figure 8. Further CO line profiles, averaged over each square degree from I = 325-330°, in units of Tmb (K)- Other details are as in 
Figure^ 
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Figure 9. intensity maps with contours in Kkms“^, corrected for the beam efficiency ryxB = 0.55, from v = —110 to 

—60kms“^ integrated over 10kms“^ intervals. The contours are at 5, 15, 25, 35 Kkms“^. Note that the boundary of the survey 

is evidenced by spots of poor signal to noise along the b = +0.5° boundary. 


lated using the dilated mask code developed for the 


MAGMA survey of the Magellanic Cloud s in CO (|Wong 
et al.|2QlT]). The dilated CPROPS mask (Rosolowsky & 


Leroy II 2006 ) is calculated by identifying regions of high 


significance (> 4cr) and expanding these to connected 


regions of lower significance (> 3cr). Masking the data 
in this way reduces the number of high noise artefacts in 
the intensity images, although some remain. In particu¬ 
lar, the upper boundary of the main survey at I = 320- 
327°, b = +0.5° is clearly visible as a series of high noise 
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Figure 10. As per Figure]^ intensity maps with contours in Kkms“^, from v = —60 to —10kms“^ integrated over 

10kms“^ intervals. Note that from the third map the intensity scale differs from Figuredue to the weaker fluxes in these velocity 
ranges, and the contours are now at 1, 8, 15, 22 Kkms“^. Due to the weaker signal fluxes there are more localised instances of 

high noise apparent in the contours in this Figure, most obviously the column at I 323.8° in the lower three maps. 


points in both the intensity map and the 

contours. 

To aid interpretation of the data, Figure shows 
derived relationships between the radial velocity and 


galactocentric radius (i?, left) and the distance from the 
Sun right) for the central longitude of each square 
degree in DRl. This is an expansion of Fig. 9 in Paper I, 
now showing the variation across the whole ten degrees. 
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Figure 11. As per Figure intensity maps with contours in Kkms from v = —10 to +40kms ^ integrated over 

10kms“^ intervals. 


As in that figure, the right plot shows the near-far am¬ 
biguity for negative velocities between 0kms“^ and the 
tangent velocity along each sightline (decreasing with 
/). From the intensity maps in Figure]^ it is clear that 
a number of molecular clouds occur at ‘forbidden’ ve¬ 
locities (possessing velocities that are not possible un- 
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der the assumption that the rotation curve holds, i.e. 
< —95kms“^); in this analysis any such velocities are 
assigned to the tangential distance. Clouds with posi¬ 
tive velocities have source distances greater than twice 
the tangent distance. 




































Mopra CO Data Release 1 


11 


I = 320.5-329.5 degrees 


I = 320.5-329.5 degrees 




Figure 12. Radial velocity-distance relationships calculated for the centre of each square degree in DRl using the |McClure-Griffiths 
fc Dicke^ (|2007| > rotation curve for the inner Galaxy (i.e. negative velocities, with R < Rq) in the fourth quadrant, with the [Brand fc 
Bhtz| ( |1993| > curve for the outer Galaxy (i.e. positive velocities, and scaled to give the same orbital velocity at R = Rq)- To the left, the 
galactocentric radius in kpc is plotted against radial velocity, Vlsr in kms“^. To the right, distance from the Sun, R, in kpc is plotted 
against Vlsr- Near-distance solutions assume D < Rtangent 7.0kpc. Far-distance solutions are for Rtangent < D < 2F)tangent- 


Position—Velocity 



Galactic Lonaitude 

Figure 13. position-velocity image for DRl, with the Galactic longitude, I on the x-axis, against the Rlsr radial velocity in 

kms“^ on the y-axis. The data have been averaged over the central degree in latitude; residuals from poor data are evident in the 
vertical striping. The solid lines are model positions of the centre of the spiral arms in a four-arm spiral galaxy with pitch angle of 
12.5°, a central bar length of 3 kpc and a Sun-Galactic centre distance of 8.0 kpc (parameters from |Vallee|2014| . From the lower edge 
on the left axis these are: the Norma arm (near at — 100kms“^ and far at —60kms“^; black dotted line), the Scutum-Grux arm (near; 
red dot-dashed line), the Sagittarius-Garina arm (near; blue dashed line), the Scutum-Crux arm (far; red dot-dashed line) and the 
Sagittarius-Carina arm (far; blue dashed line). Only a few clouds are visible in this furthest arm, and only in ^^CO and ^^CO. 


Figureshows a position-velocity slice of the 
emission averaged over the primary degree of latitude of 
the DRl data cube, overlaid with model positions of the 
spiral arms in a four-arm spiral galaxy. The parameters 
of the model galaxy have been taken from Vallee (2014) 


and depict a four-arm spiral galaxy with a pitch angle 
of 12.5°, a central bar length of 3 kpc and Sun-Galactic 
centre distance of 8.0 kpc. While the molecular clouds 
clearly lie in striped bands across the position-velocity 
diagram, the exact centre of the spiral arms is harder 


PASA (2015) 

doi:10.1017/pasa.2015.xxx 

























12 


Braiding et al 


Position—Velocity, b=0.45 




Figure 14. position-velocity images for DRl, with the 

Galactic longitude, I on the ar-axis, against the Rlsr radial ve¬ 
locity in kms“^ on the y-axis. The data in each plot have been 
averaged over 6 arcmin in latitude, the upper image is centred on 
b = 0.45°, the central plot is the midplane 6 = 0° and the lower 
image is centred on 6 = —0.45°; residuals from poor data are ev¬ 
ident in the vertical striping. 


to determine. The meta-analysis of Vallee (2014) calcu¬ 
lated the radius of the spiral arms to be near 400 pc from 
mid-arm to dust lane, which suggests that it is difficult 
to identify clearly the arm in which individual clouds re¬ 
side, particularly in regions such as I = 329-330°, where 
the Norma and Scutum-Crux arms are close in velocity 
space. 
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The distribution of molecular clouds with latitude 
is illustrated in Figure which shows the position- 
velocity distribution of emission vertically av¬ 

eraged over the 6 arcmin surrounding each of 6 = 0° 
and 6 = ±0.45°. The differences between these images 
is striking, showing how discrete the molecular gas is 
even within the Galactic plane. There is very little CO 
located at the ‘forbidden’ velocities indicated by Figure 


12 at 6 = —0.45°, while the fastest negative velocities 


occur at the higher latitudes. It is also worth noting 
that the vertical streaking, caused by high noise resid¬ 
uals in the data, occurs at different longitudes in the 
three position-velocity diagrams; these suggest that the 
methods used to clean high-noise columns work well, 
although individual pixels with poor baselines may re¬ 
main. 

Figure p±| shows the distribution of mass in DRl, un¬ 
der the assumption that the clouds with negative veloc¬ 
ities are located at the near distance (providing a lower 
limit on the true molecular mass), and that the X-factor 
used to link the CO intensity to the H 2 column den¬ 
sity is Xco = 2.7 X 10^^cm“^ (Kkms“^)“^. This value 


is adopted for |6| < 1° using figure 11 of Dame, Hart- 


|mann fc Thaddeus (2001), and while it is larger than 
their value of Xco = 1.8 x 10^^ cm“^ (Kkms“^)“^ for 
|6| > 5°, it is within the uncertainties recommended by 
[Bolatto, Wolfire fc Leroy (2013) in their meta-analysis 
of X-factor values from observational and theoretical 
studies of the Milky Way. These simplifications were 
also adopted in Paper I to provide a lower limit on the 
molecular mass in G323; however it was also shown that 
one of the five features in the spectrum was likely mis- 
identified as being at the near distance rather than the 
far distance. This leads to an underestimate of the total 
mass of about one quarter, which is within the uncer¬ 
tainty of using a constant X-factor mass calibration. 

The mass distribution in Figure while coarse, 
shows how the lower mass limit for each square de¬ 
gree varies over these 10° in longitude. The contribu¬ 
tions to the total mass per degree from clouds at posi¬ 
tive (far distances) and negative (nearer) velocities can 
be related to features in the position-velocity diagram 
in Figure In particular, the lines of sight through 
I = 324-326° possess a large component of mass at pos¬ 
itive velocities, arising from the clouds at ±35kms“^ 
from the far Sagittarius-Carina arm (see Figure 11), 
while the increase at I = 320° is due to a large cloud 
at u = [—10, ±10] kms“^, which cannot be associated 
with a single spiral arm using only the CO emission. The 
mass increases from a minimum at I r\j 322° of - 10® M© 
to a maximum at I r\j 328° of ^ 5 X 10^ Mq, with a to¬ 
tal of ^ 4 X 10^ Mq within the 10° main survey span, 
where |6| < 0.5°. Of this, approximately 15% is situated 
beyond the Solar circle. 
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Mass Distribution 



Longitude 

Figure 15. Mass distribution in DRl, calculated from 
the emission per square degree using Xqq = 2.7 x 

10^® cm“^ (Kkms“^)“^ (see text). The dotted line is the mass 
contribution from molecular clouds at negative velocities (as¬ 
sumed to be at the near distance), the dashed line that from 
positive velocities (sited beyond the Solar circle), and the solid 
line is the total (all in units of Solar mass). 


3.1 Sample Science 


The principal intention of this paper is to describe the 
CO data cubes that make up DRl. These may be used 
for a variety of applications, as they provide the distri¬ 
bution of the molecular gas along a substantive portion 
of the fourth quadrant of the Galaxy, as well as its dy¬ 
namical motions. From this, column densities, optical 
depths and the three dimensional molecular mass dis¬ 
tribution can be inferred, such as were outlined above. 
Comparison with infrared to millimetre-wave contin¬ 
uum images, for instance, would allow star formation 
efficiencies to be determined, and the variation of the 
Xco conversion factor between CO flux and H 2 column 
density to be examined on a galactic scale. Three sample 
science investigations that have been undertaken using 
the Mopra CO survey data are illustrated here. 

The first compares the CO emission in the G328 re¬ 
gion with that of [Cl] 809 GHz (from the HEAT tele¬ 
scope in Antarctica) and HI 21cm (from the South¬ 
ern Galactic Plane Survey; |McClnre-Griffiths, et al. 


2005) emission, discovering a cold, filamentary molec¬ 


ular cloud nearly 1° in length but only 5' across (Bur¬ 
ton et al. ]M4l ). The emission was confined to a narrow 
(2kms“^ FWHM) velocity range, with the hydrogen 
seen in self-absorption (i.e. HIS A), apparently envelop¬ 
ing the molecular gas. The authors hypothesised that 
the formation of a molecular cloud was being witnessed, 
with the molecular material possibly condensing out of 
the atomic substrate. Further analysis of the entire data 
cube, not just the narrow velocity range associated with 
this filament, is now being undertaken to investigate 
the variation of the [Cl] to CO ratio, and to search 
for evidence of ‘dark’ molecular gas. This is molecular 
gas where its normal tracer (CO) is weak, or absent 


(i.e. H 2 without CO), so that the molecular gas has not 
been detected, or at least the amount has been under¬ 
estimated. However if CO is absent (presumably due to 
photodissociation by far-UV photons), the carbon will 
be present as either C or C+. These species emit in the 
THz spectral regime, where the HEAT telescope oper¬ 
ates. Analysis of the complete data set (Burton et al. 
2015, in preparation) indicates that at the edges of the 
molecular clouds there is a tendency for the [Cl]/CO 
ratio to rise, as the CO itself falls, consistent with the 
existence of dark molecular gas. 

The second example relates the Mopra CO data to 
X-ray emission. A bright X-ray fiare was observed from 
the binary star Circinus X-1 in late 2013 by three space¬ 
craft {Swift, Chandra and XMM-Newton). Over a three 
month period four well-defined rings were seen in the X- 
rays, extending outwards with time to a distance of 14' 
from the central source. This is in the G322 portion of 
our CO survey. Comparison of the Mopra data with the 
images from the X-ray satellites showed that each of the 
X-ray peaks in the rings is correlated with correspond¬ 
ing emission in the molecular gas. However, each has 
a different velocity, ranging from —81 to —32kms“^ 
(see [Heinz et al. 2015). The X-rays are being scat¬ 
tered off clumps of dust associated with the molecular 
gas to produce the light echoes. Using the distances to 
these clumps inferred from the Galactic rotation curve, 
and the scattering geometry a direct kinematic distance 
could be determined for the source. This was found to 
be 9.4 kpc, significantly greater than the previously esti¬ 
mated 4kpc. It implies that Circinus X-1 is a frequent 
super-Eddington source, producing a narrow jet with 
a collimation angle of 0 < 3° and a Lorentz beaming 
factor of T > 22. 

The final example relates the Mopra CO data to an 
unidentified radio continuum source known as PMN 
1452-5910, thought possibly to be related to an SNR 


evident at 843MHz from MOST data (Bock, Large & 
Sadler|[T9^ , as analysed by [Jones fc Braidingj ( |2Q1^ 


Three features are apparent in the Mopra CO spectrum 
towards this source, at —50, —42 and —lkms“^. Their 
morphology, however, in comparison to that of the ra¬ 
dio continuum, suggests that it is the —lkms“^ CO 
feature that is associated with the radio, so providing 
a distance to the source. This also identifies the source 
with an H 2 O maser at a similar velocity, and enables an 
interpretation of infrared images from the Spitzer tele¬ 
scope to be made. Jones & Braiding concluded that the 
source is a previously-unidentified massive star forming 
region about 13 kpc away. They then determine param¬ 
eters for the region: ^2x10^ Mq of molecular gas with 
an average column of ^ 10^^cm“^, enclosing several 
UCHll regions with emission measures ^ 10^“^ pc cm“^ 
and electron densities r\j 10^ cm Three OB stars (2 
of type BO, 1 of 08) can account for the radio emission. 
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4 CONCLUSIONS 


In their forthcoming Annual Review of molecular clouds 


in the Milky Way, Heyer & Dame (2015) state: “For fur¬ 


ther progress in [molecular cloud census studies], sur¬ 


veys comparable to the Galactic Ring Survey ( Jack- 
son et al. 2QQ6[ ) in resolution and sensitivity in the 


fourth quadrant are required.” The first 10° of such a 
survey have been presented here. The Cherenkov Tele¬ 
scope Array will produce a new survey of the high en¬ 


ergy gamma-ray sources in the Galactic plane (Dubus 


et al. 2013) at a resolution approaching 1 arcminute 
(Bernlohr et al. 2013). High resolution molecular gas 


surveys such as this will be required to properly identify 
these sources amongst the diffuse TeV gamma-ray emis¬ 


sion components CTA will likely detect (e.g. Acero et al. 


2013). The Mopra CO survey of the southern Galactic 


plane aims to address these community needs by observ¬ 
ing the entirety of the fourth quadrant along the plane 
from I = 270-360° where |6| < 0.5°, at a spatial resolu¬ 
tion of 35 arcsec and a spectral resolution of 0.1 kms“^, 
in the J = 1-0 lines of ^^CO, ^^CO, C^^O and C^^O. 
This is a significant improvement over previous studies 
conducted of the fourth quadrant of the Galaxy. 

Data Release 1 from the Mopra CO survey covers 
the 11.5 square degrees from I = 320-330°, b = ±0.5° 
and I = 327-330°, b = ±0.5-1.0°. Included in this pa¬ 
per are a number of metrics describing the quality of 
the dataset, which has a typical rms noise of ^ 1.3 K in 
the ^^CO frequency band and ^ 0.5 K in the other iso- 
topologue frequency bands, and a discussion of the im¬ 
provement in the data taken during the winter months 
compared to the initial square degree observed in sum¬ 
mer that was published in Paper I. The line intensities 
in the Mopra CO survey are ^1.4 times those of the 


Dame, Hartmann & Thaddeus (2001) survey, and our 


higher resolution CO line profiles match those of their 
survey quite well. Sample line profiles for each square 
degree are presented, as well as integrated intensity 
maps that show the filamentary nature and complex¬ 
ity of the molecular clouds in the survey. The primary 
10° X 1° strip of the Galactic plane in DRl has been 
found to contain ~ 4 x 10^ Mq of molecular gas. 

The Mopra CO survey is ongoing, with the intention 
of covering the entirety of the fourth quadrant of the 
Galactic plane by the end of 2015. The data are be¬ 
ing made available at the survey websit^ and in the 
CSIRO-ATNF data archiv^ as they are published. 
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